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The effect of an external force on the kinetics of diffusion-assisted escaping of Brownian particles 
from a potential well is analyzed in detail. The analysis is made within the two-state model of 
the process which is known to be valid in the deep well limit in the absence of external force. 
The generalized variant of this model, taking into account the effect of the force, is shown to be 
quite accurate as well for some shapes of the well. Within the generalized two-state model simple 
expressions for the well depopulation kinetics and, in particular, the for the escape rate are obtained. 
These expressions show that the effect of the force (F) manifests itself in the escape rate dependence 
on the only parameter ip = Fa/(2kbT), where a is the Onzager radius of the attractive part of the 
well U(r), defined by the relation \U(a)\ = kbT. The limiting behavior of this dependence in the 
cases of weak and strong force is analyzed in detail. Possible applications as well as the relation of 
the results of the analysis to those obtained earlier are briefly discussed. 
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I. INTRODUCTION 

The effect of external force on mechanisms and kinetic 
properties of condensed phase diffusion-assisted reaction 
processes is considered in a large number of works both 
experimentally and theoreticallyfA±S The active interest 
of scientists to this phenomenon results from its great 
practical importance. 

One of the most important systems, in which the ex- 
ternal force effect is investigated very thoroughly, is re- 
combining geminate ion pairs, undergoing relative diffu- 
sion in the external electric field^^ Most of theoretical 
studies analyze the kinetics of the processes within the 
simplest model, which reduces the problem to solving 
the Smoluchowski equation for probability distribution 
function (PDF) of particles diffusing in a pure Coulomb 
potential (with an external force) and reacting with the 
rate highly localized at short distances!^ Even in this 
most simple formulation the problem can, in general, 
be solved only numerically, though detailed analytical 
analysis of some simple variant of the problem have also 
been made^^ for example, within the prescribed diffu- 
sion approximation^ 

Recent advances in time resolved investigations of 
charge transfer and escaping processes in fast geminate 
reactions- and, in particular, geminate recombination of 
ion pairs in polar liquids^&^ i 10 ' 11 ! 12 ' 13 inspire further de- 
velopment of theoretical methods of the analysis of the 
considered problem. The main challenge of the theo- 
retical studies consists in the correct description of the 
manifestation of specific features of the interparticle in- 
teraction (in real liquids) in the reaction kinetics in a 
tractable form simple enough to be suitable for applica- 
tions. 

In the majority of above mentioned theoretical works 
no specific features of the form of the interaction poten- 
tial for the probe (Brownian) particles at short distances 
(of order of molecular size) have been taken into account. 
In the condensed phase, however, the distance depen- 



dence of the potential at short interparticle distances r 
can be strongly modified by interaction of particles under 
study with those of the medium^ This modified interac- 
tion is usually characterized by the so called mean force 
potential (MFP), which in a physically reasonable form 
incorporates the medium effect and, in particular, dis- 
creteness of the medium at short distances. The interac- 
tion with the medium particles is known to result in the 
wavy behavior of the MFP at short distances. Moreover, 
in some systems the medium effect results in the well- 
type shape of the MFP at short distances (see Fig. 1) 
with the a markedly high barrier at distances r of order 
the distance of closest approach d. This effect is found, 
for example, in the case ion pairs in polar liquids i^ill 

Concerning the applicability of well-type approxima- 
tion for the real MFP shape, it is also worth mentioning 
the additional reasoning: from mathematical and kinetic 
points of view any attractive potential can be considered 
as well-shaped in the absence of (or low) reactivity of par- 
ticles at r ~ d. The only difference of this type of wells 
from those shown in Fig. 1 is in their urge-like shape at 
r ~ d. 

The well-type shape of the MFP (with the reaction bar- 
rier at r ~ d) results in the formation of the quasiequilib- 
rium state within the well, which can also be considered 
as a cage state. In the absence of external force the ki- 
netics of diffusion-assisted depopulation of the initially 
populated cage state is analyzed in detail in a number of 
papers i 15 i 16 ' 17 i 18 In the limit of deep well depth the prob- 
lem is shown to be accurately described with the model 
of two kinetically coupled states: the quasiequilibrium lo- 
calized state within the well and the free diffusion state 
outside the well>il 

The two-state model enables one to obtain the well de- 
population kinetics in a relatively simple analytical form. 
This kinetics, determined by the monomolecular reactive 
passing over the barrier at r ~ d and escaping from well 
(cage), appears to be non-exponential, in general. In the 
limit of deep well, however, the deviation from the expo- 
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FIG. 1: The picture of the interaction potential u(r) for two 
models of its shape: narrow (dashed) and wide (full) well; 
a„ and CL-w £1X6 the Onzager radii for these models [defined by 
u(a v ) = 1, (y = n,w)], d is the distances of closes approach, 
and rb is the coordinate of the bottom. 



of which models the reaction within the well. 

The main purpose of the work is to analyze the ki- 
netics of diffusion-assisted escaping from the well in the 
presence of the external force — F, i.e. escaping from 
the well of the potential Uf = U(r) + (F • r). For 
definiteness the force is assumed to be directed along 
the axis z: F = (0,0, F). The analysis can conve- 
niently be made in spherical coordinates in which r = 
(r sin 9 cos </>, r sin 9 sin </>, r cos 9) . 

The diffusive space-time evolution of the Brownian 
particle is described by the PDF p(r, t) = p(r, 9, (f>\t). In 
general, in spherical coordinates the PDF depends on all 
three variables r, 9, and 0. However, in the considered 
case of isotropic diffusion and isotropic initial condition 
p(r, t — 0) = Pi(r) the PDF p(r,t) is independent of the 
azimuthal angle <fi so that p(r,t) = p(r, 9\t). In our anal- 
ysis we will assume that particles are created within the 
well at initial distance ~ r^: 



nential kinetics is shown to be fairly smalh 16 i 

In this paper we generalize the two-state model to de- 
scribe the effect of the external force on the well depopu- 
lation kinetics. The effect is assumed to result only from 
the force induced change of the rate of escaping from the 
well and will concentrate on the discussion of the mani- 
festation of the force in the escaping kinetics. 

Within the generalized two-state model we derive sim- 
ple formulas for the well depopulation (or escaping) ki- 
netics for different models of the well shape. Moreover, 
taking into account that in the limit of deep well (which is 
of main interest of the discussion) the kinetics is close to 
exponential, in our study we restrict ourself to the anal- 
ysis of the escape rate, only. The analysis shows that the 
escape rate is strongly affected by the external force and 
the specific manifestations of this effect depends on the 
well shape. The force effect can be characterized by the 
rate dependence on the only parameter. In the cases of 
weak and strong force the limiting analytical expressions 
for this dependence are obtained and briefly discussed. 



II. FORMULATION OF THE PROBLEM 

We consider the dynamics of the Brownian particle un- 
dergoing diffusive motion in a three dimensional spher- 
ically symmetric MFP well U(v) = U(r) centered at 
r = 0, where r is the vector of position of the parti- 
cle. The dependence of U(r) on the distance r = |r| is 
schematically shown in Fig. I. The MFP will be char- 
acterized by three parameters: the distance a (called 
hereafter the Onzager radius), at which U(a) = k B T 
[in Fig. 1 a n and a w denote the Onzager radii, corre- 
sponding to two models of the well shape, discussed be- 
low (in Sec. V)], the distance of closest approach d, and 
the radius of the bottom of the well whose energy is 
U{r b ) = -U b = -(k B T)u b . At r ~ d the MFP U(r) is 
assumed to be of the shape of a barrier, diffusive crossing 



p(r,t = 0) = (47rrl)- 1 S(r~ n ), (2.1) 

The PDF p(r,9\t) satisfies the Smoluchowski equation 

p = V r [D(r)(V T p + pVrUf)}, (2.2) 

where V r is the gradient operator, 

u,(r) = u(r) + (fT), (2.3) 

with u(r) = U(r)/(k B T) and f = F/(k B T), is the di- 
mensionless MFP, and D(r) is the diffusion coefficient, 
which in our analysis is assumed to be independent of 
r: D(r) = D [though some possible effects of D(r)- 
dependence can also be studied]. 

Note that the reaction kinetics for pairs of interacting 
Brownian particles, say a and b, is described by the equa- 
tion similar to eq. (|2.2p with r = r a — and parameters 
expressed in terms of those for separate particles. 2 > 4 

In the absence of force the kinetics of escaping from the 
spherically symmetric short range potential well u(r) is 
analyzed earlieri 15 i 16 i 17 i 18 Here we extend the approach 
applied in these works to describe the effect of external 
force. This approach is based on the approximate solu- 
tion of the eq. (|2.2|) in the limit of deep well, in which the 
solution can be obtained in analytical form by expansion 
in a small parameter T r /r e I, where r r ~ (a — d) 2 /D is 
the time of equilibration within the well and r e ~ r r e _tl ° 
is the time of escaping from the welh 16 ' 17 

Analysis of this solution showsi 7 - that in the lowest 
order in the parameter r r /r e the Smoluchowski approxi- 
mation (|2.2[) is equivalent to the model of two kinetically 
coupled states: the state within the well and the state of 
free diffusion outside the well. To clarify the main points 
of this model we will first briefly consider the case of the 
absence of the force F = 0. 
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III. TWO-STATE MODEL IN THE ABSENCE 
OF FORCE 

Originally, the two-state model was proposed to treat 
the kinetics of the diffusion-assiisted escaping from the 
well in the case F = 0, in which the PDF p(r, t) depends 
only the distance r — |r|: p(r,t) = p(r,t). 

Analysis show s 16 ' 17 that in this case in the lowest or- 
der in r r /r e <C 1 the escaping kinetics can be described 
within the model of two kinetically coupled states: highly 
localized state within the well and free diffusion state out- 
side the well. The evolution of the state within the well 
(d < r < a) is determined by the well population 



i(t) = 4tt / 

Jd 



dr r 2 p(r, t), 



(3.1) 



while the evolution of the state outside the well (r > a) 
is governed by the distribution function c(r, t). The func- 
tions n(t) and c(r,i) satisfy simple kinetic equation s 16 ' 17 



h = S t K + c(a) - (K- + w r )n (3.2a) 
c = L r c+(Sf x K-n-K + c)8(r-a), (3.2b) 

which should be solved with the initial condition 



n(0) = l and c(r, 0) = 



[implied by eq. (|2.ip ] and the boundary conditions for 
c(r,t) written as V r c(r, t)\ r=a = (this condition corre- 
sponds to reflection at r = a) and c(r — * oo) = 0. . 
In eqs. (|3~2"j) Si = 4ira 2 and 



L r = Dr~ 2 V r (r 2 V r ) 



(3.4) 



is the radial part of the free diffusion operator, in which 
V r = d/dr. The terms proportional to K± describe the 
above-mentioned kinetic coupling (transitions) between 
the state within the well, located at r = a, and the free 
diffusion state outside the well. In the considered limit 
r r /r e <C 1 the transition rates K± satisfy the relations: 17 

K ± ^oo and K_ / K + = K° e = a 2 / Z Wl (3.5) 

where 



drr 2 e- u{r) 



(3.6) 



d<r<a 



is the partition function for the well u(r) — f7(r)/(fcsT). 
Formula (|3.5[) . obtained by comparing equilibrium solu- 
tions of eqs. (|2.2[) and §j£2§ r^ represents the detailed 
balance relation for transitions between the states of the 
two-state model. 

The term w r n(e) in eq. (|3.2ap describes the effect of 
the first order reaction (in the well) with the rate 



= [D/Z w ) 



drr~ 2 e u{r) 



rr^d 



(3.7) 



Solution of eqs. (O yields^ 6 -! 7 - 

■ioo+0 



nit) = ± J de 



exp(et) 



— 2OO + 



e + w r + a 2 K e V{e) 



(3.8) 



The function V(e) is directly related to the Green's func- 
tion of the operator controlling diffusion outside the well 
[with the reflective boundary condition V r g\ r=a = 0] 16 i 17 



G(r,n\e) = (r\(e-L r )- 1 \r l ) 



(3.9) 



V{e) = 1/G(a, a\e) = D[ a - 1 + (e/L>) 1/2 ]. (3.10) 



Substitution of formula (|3.10p into eq. (|3.8p leads to 
the following formula for the well population n(t): 17 



n(t) 



1 



100+0 d£ exp[e(w t)} 



— 2OO + 



1 + e + 7£ 



1/2 



where 7 = {w e /wQ) 1 / 2 (a 2 w e / D) 1 / 2 . The rate 



w = w e + w r 



(3.11) 



(3.12) 



(3.3) is a sum of the rate of escaping from the well 



= f~ (/ ° drr ~ 2eM(r) ) = ° a/Zw (3 ' 13) 



and the rate of reaction in the well w r [see eq. (|3.7[) ]. 

The function n(t) (|3. 1 1|) is, evidently, non-exponential. 
Its analytical properties are essentially determined by 
the parameter 7 w a/In, where Id = (D/wq) 1 ^ 2 is the 
average distance of diffusive motion during the lifetime 
To = Wq 1 of the particle in the well. It is easily seen 
that for deep wells Id ~> a, and therefore 7 < 1. The 
parameter 7 controls the onset of the change of the func- 
tional form of n(t): n(t) = exjp(—wot) air < 111(1/7) an d 
n(t) - l/t 3 / 2 at r > 111(1/7). 

It is worth emphasizing that despite the complexity of 
the escaping kinetics n(t), in general, it is quite close to 
the exponential for deep wells (for 7 <ti \), as it follows 
from the above expressions. Moreover, independently of 
the well depth average steady state characteristics of type 
of the time averaged kinetic parameters just coincide with 
those corresponding to the exponential part of the kinet- 
ics. For example, the lifetime of particles in the well 



To 



dtn{t) 



(3.14) 



This means, in particular, that the steady state char- 
acteristics are quite informative and can be used in the 
discussion of the escaping kinetics in the more compli- 
cated case of the presence of an external forces. 
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IV. TWO-STATE MODEL IN THE PRESENCE 
OF FORCE 

A. Kinetic equations 

In the presence of force [f = F/(fcsT) ^ 0] the poten- 
tial Uf(r) = u(r) + (f • r) in eq. (12. 2p . is not isotropic, 
which results in the dependence of the solution p(r, t) on 
particle polar angle 9: p(r,t) = p(r,6\t). 

It is important to note that in the case f ^ the two- 
state model is also valid for a variety of shapes of the 
potential well u(r) though some additional analysis of 
the corresponding validity criteria are certainly required. 

Similarly to the case of the absence of force, in any 
variant of the anisotropic two-state model (for f ^ 0) 
the kinetics is also described by two functions: the pop- 
ulation n(8\t) within the well and the PDF c(r,6\t) of 
particles outside the well. These functions, however, de- 
pend on the polar angle 9. Within the two-state model 
the effect of the force f shows itself in the dependence of 
kinetic parameters on 9. The form of this dependence is 
determined by the particular variant of the model. 

In general, two-state kinetic equations, describing evo- 
lution of PDFs n(9\t) and c(r,9\t) in the presence of an 
external force, can be written by analogy with eqs. (|3.2p : 

h = SiK+(9)c(a\t) + [L c - (K_(9) + w r )]n, (4.1a) 
c = L f c+[S- 1 K^{9)n- K + (9)c]S(r - a), (4.1b) 

where 

Lf = DVr{Vr + f) (4.2) 

is the operator, describing diffusion outside the well, and 
L c is the Smoluchowski operator in {#}-space which con- 
trols orientational relaxation of the PDF in the well. 

The essential difference of anisotropic equations from 
isotropic ones consists in the orientation dependence of 
rates, K + (9) and K-{9). Below, in accordance with eq. 
p.5p . in the considered limit r r /r e <C 1 we will assume 
the transition rates K± to satisfy the relations^ 

K ± -> oo and K-(9)/K+(6) = K e {9). (4.3) 

Therefore in this limit ^-dependence of rates K± (9) show 
itself in that of the the equilibrium constant: K e (9). The 
form of the function K e {9) is determined by the shape 
of the well. Some model well shapes and corresponding 
K e (9) dependences, as well as applicability of the corre- 
sponding two-state models, are discussed below. 

Equations (|4.1[) should be solved with boundary con- 
ditions 

(V r + f cos 9)c\ r=a — and c\ r ^oo — > 0, (4.4) 

first of which describes reflection of particles (diffusing in 
the state outside the well) at r = a. The initial condition 
is assumed to be isotropic and given by eq. (|2.ip . 



In what follows it will be convenient to represent func- 
tions n{9\t) and c(r, 9\t) in the form of vectors |n(t)) and 
|c(r, £)), whose components are obtained by expansion 
of these functions in the orthonornmal basis of properly 
normalized Legendre polynomials (spherical functions) 19 

= (Z + l/2)P;(cos0), (4.5a) 

(Yi\ = [ d9 sin 9Pi{cos 9)... (4.5b) 
Jo 

with (1 = 0,1,...): 

OO OO 

\n)=Y t m\Yi) and |c)=^q|^), (4.6) 

(=0 1=0 

where for any vector |x(#)), (x = n, c), its components 
Xi are defined by 

Xl = {Yi\x)= I d9 sin 9 P l (cos9) X (9). (4.7) 
Jo 

In term of this vector representation the initial condi- 
tion can conveniently be written in the form, explicitly 
displaying its independence of orientation: 

\p l )^(2nrf)- 1 \Y )S(r-r l ). (4.8) 

As for the initial condition, it is worth noting, in ad- 
dition, that in the most realistic limit of orientational 
relaxation within the well much faster than the escaping 
from the well the escaping kinetics is insensitive to the 
orientational dependence of the initial condition. 

B. Escaping kinetics 

Equations (|4.ip can be solved by the method, applied 
above in the case of the absence of force, but with the use 
of expansion of n(9\t) and c(r, 9\t) in spherical functions 
\Yi), i.e. vector representation \n(t)) and |c(r|t)) [see eq. 
(I4.6p ], The solution yields for the Laplace transform 

POO 

|5(e)) = / dte-*\n(t)) : (4.9) 
Jo 

|n(e)) = (e + w r + L c + a 2 V{e)k e y 1 \n l ), (4.10) 

where for the initial condition (|4.8p |n$) = (l/27r)|lo). In 
this expression the equilibrium constant K e is the oper- 
ator, which indicates its dependence on the orientation 
angle 9, and 

V = G- 1 (a,a\e) (4.11) 

is the operator in the space of |Yi)-states, expressed in 
terms of the evolution operator for diffusive motion out- 
side the well (evaluated at r = r{ = a): 

G(a,a\e) = (a\(e — I//) _1 |a) 

= e- vcos6 {a\(e- A / )- 1 |a)e vcose , (4.12) 
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in which tp = fa/2 and 

\ f = D(L r + r- 2 L e -y 2 ), (4.13) 

with 

L fl = (l/Bto0)Vfl[sin0(V fl )]. (4.14) 

In the expression (|4.14[) L r is the operator of radial 
diffusion, defined in eq. (|3.4p , and Lg is the operator of 
free orientational diffusion, diagonal in |Yj)-basis: 

Le = -YZo l{l + 1)lYl){Yl1 (4 - 15) 

In what follows we will mainly restrict ourselves to the 
analysis of the total population of the well no(i), for the 
Laplace transform of which one gets the expression 

n (e) = 2tt / d0 sm6n(6, e) = 27r(F |n(e)) 
Jo 

- (K |[e + uv+£c + a 2 y(e)£ e ]~Vo>> (4.16) 

Moreover, since the escaping kinetics appears to be 
fairly close to the exponential (see Sec. Ill), to charac- 
terize the process it is sufficient to calculate the mean 
inverse lifetime tq of particles in the well defined by eq. 
(pU4l) : 

f =n (0) = (Y \[w r + L c + a 2 V(0)k e ]~ 1 \Y ). (4.17) 

C. General formulas 

Formulas (|4.12I(|4.14p allow us to evaluate the opera- 
tor V(e) in analytical form and, therefore, analyze the 
behavior of no(e) and fo relatively easily. 

In the evaluation it is worth taking into account the 
useful relation which simplifies the differential operator 
in the radial space: 

(a\(e-A f )- 1 \a) = {a\(e-\ f )- 1 \a), (4.18) 

where 

\ f = D{V 2 r + r- 2 L e -\f). (4.19) 

The evolution operator (a|(e— A/) _1 |a) can be obtained 
in analytical form 2 ^ with the use of two linearly indepen- 
dent operator solutions ip-(r) and ip+(r) of equation 

(e - h)i>± = (4.20) 

in which the operator Lg is treated as a parameter. These 
solutions satisfy two boundary conditions corresponding 



to those given in eq. (|4.4[) [after change of variable c(r) = 

e -(/rcoBfl)/2^ r )] 

(V r + i/w)^_| r=0 = and ^+1^00 -» 0, (4.21) 
where 

oo 

u> = \Yi)(Yi\ cos 6\Y m )(Y m \ (4.22) 

i,m=0 

is the matrix representation of the function cos 9. The 
matrix elements (V; | cos 6\Y m ) = (Yi\Pi(cos8)\Y m ) are 
evaluated analytically 19 though the corresponding for- 
mulas will not be needed in our further analysis. 

Both solutions ip+ (r) and ip- (r) are expressed in terms 
of matrices of modified Bessel functions 21 

K(r) = VtJ2Zo lYl)Kl +i ikrm > (4 - 23) 

_ x -OO 

I(r) = V^E ;=0 l y ')^+^ fcr )^l : ( 4 - 24 ) 

^+(r) = K(r), */L(r) = J(r) + K(r)fi, (4.25) 

where 

fc(e) - v/(//2) 2 + e/^ (4.26) 

and k is the matrix determined by the boundary condi- 
tion at r = a [see cq. (|4.21| : 

k = [V r £(r) - ^(r)]- x [g7(r) - V r /(r)]| r=a , (4.27) 
in which 

g = a _1 (l - ipij) = a _1 (l - </jcos0), (4.28) 

with = fa/2. 

It is worth noting that the matrices K and I do not 
commute with Co and, therefore, the order of matrices in 
the products of the matrices in expressions (|4.25|) - (|4.27|) 
is important. As a result of these special commutation 
properties of the matrices, the matrix solutions ip+(r) 
and ip-(r) do not commute either. 

The general representation of the evolution operator 
(r|(e — A/) -1 |ri) in terms of non-commuting solutions 
i/>+ (r) and i\>- (r) is proposed and thoroughly discussed 
in ref. [20]. This representation generalizes the one well 
known in the case of scalar solutions ip+(r) and ip-i r )- 
In general, the proposed representation is fairly compli- 
cated, however, in the particular case of solutions given 
by eqs. (|4.23j) and (I4.24| it reduces to a more simple one: 
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(r|(e - A/) _1 l r i) = -D-^KtfliriWair - n) + i{r)K{ri)9 H (ri - r) + K(r)kk(r t )]/W{Kj). (4.29) 



In this expression 6h(x) is the Heaviside step function 
and 

W(K, I) = V r K(r)I(r) - VJ{r)K{r) (4.30) 

is the Wronskian of two solutions which is a scalar func- 
tion [W(K,I) = — 1]. The validity of the expression 
(|4.29p can easily be verified by direct substitution of it 
to equation similar to eq. (|4.20[) but with delta-function 
in the right hand side. 

For the particular case r — ri — a formula (|4.29|) yields 

(a\(e-\ f )- 1 \a)=D- 1 [q-q K (e)}-\ (4.31) 

where 

oo 

q K (e) = V r k{r)/K{r)\ r=a = £ \Y x )q Kl (e)(Y l \ (4.32) 

1=0 

with 

q Kl (e) = {I I a) + fc(e)*Q_ i (k(e)a)/K l+ 1 (k(e)a). (4.33) 

Substituting the expression (|OT|) into eq. |4~T8)) and 
then into eqs. (|4.13p and (|4.11|) we obtain fairly simple 
formula for V(e): 

V(e) = De- vcose (q-q K )e vcose 

= D[a -1 (l -tpcos6) - e- ¥,cose q K e ¥,cose ].(4.34) 

For out further analysis of the escaping kinetics no(t) 
we need to specify of the operator L c describing orien- 
tational relaxation in the well. Naturally it should be of 
the Smoluchowski-like form: 

L c = Dcismey^elsine^e + V u)], (4.35) 

where D c is the orientational diffusion coefficient u(Q) 
is the effective orientational which is determined by the 
shape of the well (see below). 

For the sake generality, we will not assume any relation 
of the value of D c with that of the diffusion coefficient D 
outside the well (of type of D c ~ D/r%). 

Moreover, in the considered limit of large well depth it 
is quite natural to assume that orientational relaxation 
is much faster than well depopulation. 

D. Fast orientational relaxation in the well 

The fast orientational relaxation limit implies that 
D c 3> wq, where wq is the rate of depopulation of the 
well defined in eq. (|3 . 14|) . This relation means that af- 
ter some time ~ r c = D~ l of orientational relaxation (of 



I 

the initial population in the well) the vector of well pop- 
ulation |n(t) remains close to the equilibrium one \^ e ) 
during the process: 

|n(t)> « no(t)|* e >, (4.36) 

where 

|* e ) = Zg 1 e- UbW , Z g = f d9sm6e- abW . (4.37) 

Jo 

Note that within bra-ket notation the adjoint vector (ip e \ 
coincides with (Yq \ and is given by formula 

(* e | - (Y \ = f dd sin0..., (4.38) 
Jo 

which can be confirmed by the relation (^ e \L c — di- 
rectly following from the definition of L c [see eq. (I4.35|) ] . 
With the use of this formula one can easily find that \ifj e ) 
satisfies the normalization condition (\f r e |'\E r e ) = 1. 

In what follows we will restrict ourselves to the analysis 
of the escaping kinetics just in this limit of fast orienta- 
tional relaxation. 

For fast orientational relaxation the splitting SL C of 
eigenvalues of the operator L c (SL C ~ D c ) is much larger 
than o 2 !!^^!! ~ w e . In such a case in the lowest order in 
the parameter w e / D c C 1 we can significantly simplify 
the general expressions (14.16[) and (|4.17|) for no(e) and the 
inverse average lifetime wo = f^ 1 = 1 (e), respectively, 
as follows: 

no(e) = [e + w r + a 2 e \V(e)K e \^ e )]-\ (4.39) 

and 

w = f- 1 = w r + a 2 (y e \V(0)k e \y e ), (4.40) 

Equation (|4. 39[) presents the main result of the work 
for the kinetics of the well depopulation in the limit of 
fast orientational relaxation. 

Formula (|4.39[) shows that, in general, in the fast ori- 
entational relaxation limit the escaping kinetics no(t) 
is non-exponential, however, as it has been mentioned 
above in Sec. Ill in the considered case of deep wells its 
deviation from exponential is very small and, therefore, 
the kinetics can quite accurately be characterized by the 
only parameter, the mean lifetime f ~ l or the correspond- 
ing mean rate wq — f^ 1 . For this reason, in what follows 
we will mainly discuss these two parameters. 

Moreover, below we will concentrate on the second 
term in square brackets in eq. (|4.40[) : 

w e = a 2 (^ e \V(0)k e \^ e ), (4.41) 
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which can be interpreted as the mean escape rate. 

To apply formulas (|4.39p - (|4.4ip one need to specify the 
shape of the potential well which determines the orienta- 
tional potential u(9) in the Smoluchowski-type operator 
L c [see eq. (|4.35|) ] and, therefore, the equilibrium state 
l^e). In our work we will consider two realistic models 
of the well shape in which relatively simple analytical ex- 
pressions for the mean rate wq = f^ 1 and, in particular, 
for w e can be obtained. 



V. APPLICATION OF GENERAL RESULTS 

In this section we will analyze the specific features of 
mean escape rate w e for two variants of the well shape: 

1) Narrow well shape (shown in Fig. 1 by dashed line 
with a = a n ), for which a — d <C d and the time of 
equilibration within the well r r ~ (a — d) 2 / D <C a 2 /D <C 
r e , where r e ~ f _1 [see eq. (|3.14[) ] is the time of escaping 
from the well; 

2) The wide well shape (full line in Fig. 1 with a = 
a w ), corresponding to a small distance of closest approach 
d -C a (in which of the main interest is the region /del 
while fa > 1). This shape is schematically shown in Fig. 
1 by full line with a = a w . 

The analysis will be made in the above-discussed limit 
of fast orientational relaxation in the well with the use 
of eq. (|4.39[) for the well depopulation kinetics. In our 
study, in addition to this formula and other ones derived 
above, we will also use the representation for e ±vcose m 
terms of expansion in spherical functions Pi(cos9)£^ 

e ± vcos e = ^giijip+i)/ I+ i(^(cosfl). (5.1) 

V iD L — ' 2 

v r 2=0 



A. Narrow- well shape 

In the case of narrow well, when a — d < d, the well is 
of the shape of attractive well layer near the distance of 
closest approach d. In this limit within the wide region 
force strengths / < 1/ (a — d) we can neglect the effect 
of the force on the radial shape of the well and take into 
consideration only the dependence of well depth U(9) on 
the orientation angle 9: 

Ub(9) pa Uf(rb, 9) ps ui, + facos9 (5.2) 

with / = |f| > 0, and the force effect on free diffusion 
in the state outside the well. In eq. (|5.ip we took into 
account the smallness of the width of the well, a — d <C d, 
which leads to the high accuracy of the relation fr^ pa fa. 

It is important to note that the small value of the 
well width and, therefore, fast equilibration of the well 
population in radial direction, ensures the validity of 
the description of the kinetics in terms of the angular 
coordinate dependent well population n(9\t) introduced 



above. Noteworthy is also that the negligible force af- 
fected change of the well shape results in the absence 
of the dependence of the detailed balance relation and 
the equilibrium constant [see eq. (|3.5[) ] on the angle 
9, i.e K e = K® . In other words, the effective parti- 
tion function Z Wl which governs the constant K e , is still 
given by eq. (|3.6p . i.e. is controlled by the shape of 
the potential u(r) without external force, despite possi- 
ble strong force effect on the energy of the bottom pre- 
dicted by eq. (|5.2p . This is because the external force 
leads to the identical change of both the bottom energy 
Ub{9) pa Ufc + fa cos 9 and the energy of the free diffusion 
state at r = a: Uf(a,9) pa fa cos 9. 

The potential Ub{9) determines the kinetics of orienta- 
tional relaxation of the population in the well, which is 
described by the Smoluchowski operator (|4.35p with 

u(9) = Ub(9) — Ufc = 2(pcos9, where (p = fa/2. (5.3) 

In this case the equilibrium state within the well is writ- 
ten as 

|* e ) = Z^^)e-^ cose (5.4) 

with 

Z (cp) =smh(2(p)/cp, (5.5) 

Substitution of formulas (|5.4p and (|5.ip into the ex- 
pression (|4.4ip yields for the mean escape rate w e ((p), 
expressed in terms of the universal function Qni^)- 

w e ((p)/w° en =Q n {ip) = \ + ^coth(2^) + S n (tp). (5.6) 

In this formula ro^ = w° e is the escape rate in the absence 
of a force [see eq. (|3.13[) ]. if = fa/2, and 

oo 

S n ( V ) = 2Z 1 ^)(n/^Y / (l+ 1 2K^Mv) (5.7) 

1=0 

with Zo((fi) defined in eq. (15. 4p and 

qi(ip) = aq Kl = I + <pK t _ i ((p)/K l+ 1 (ip). (5.8) 

Formula (|5.6p shows that the force effect on the rate 
is characterized by the only parameter ip = fa/2. The 
numerical calculated universal function Q n (<p) which de- 
scribes this effect is displayed in Fig. 2a. In addition, 
some limiting specific features of the behavior of Q n (f) 
can be revealed with simple analytical expressions. 

a. Weak force limit. In the limit of weak external 
force, when p = /a/2 < 1, in two lowest orders in ip 
the mean escape rate w e can be estimated using only the 
first term (with I = 0) in the sum S n (ip) in the expression 
(EH) for Q{ip) =w e (ip)/w° n : 

Q n (<p)=QZ(<p)nl + ip = l + fa/2. (5.9) 

b. Strong force limit. In the opposite case ip = fa/2 ^> 
1, denoted as the strong external force limit, the analysis 
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2 3 4 5 <p 



FIG. 2: The force (<p = fa/2) dependence of dimensionless 
escape rates Q v {tp) = w e (tp)/w° l/ ((p) (y — n,w) [for narrow 
(n) and wide (w) wells], calculated with exact eqs. (|5.6|l and 
(|5.15p (full lines) and interpolation formulas (|5.1ip and (|5.20p 
(triangles) . 



of the dependence w e (ip) with the use of eq. (|5.7p is 
more complicated. However, the corresponding limiting 
dependence can easily be derived taking into account that 
in the limit ip 3> 1 the escaping process is, actually, one- 
dimensional, for which the <3(</?)-factor is given byiS. 



= Q S M = fa. 



(5.10) 



In deriving eq. (|5.10p we used the expression for one- 



dimensional escape rat' 



,18 



Wl 



Df/Zi, in which Z\ 



Z w /a is the one-dimensional parti- 



Jd<r<a f 
tion function. 

This relation can also be illustrated by the semiquan- 
titative estimation with the use of eq. (|5.7p by trun- 
cating the sum at I = I* ~ ^/p ^> 1. The fact is 
that the terms of the sum with I < I* give the main 
contribution to the sum since the terms with I > I* 
rapidly decrease with the increase of I. At 1 < / < I* 
we can use the estimations^ I l _i{p) ~ e v /y/(j) and 
K l _i(ip)/K l+ i(ip) w 1 so that qi<i*{<p) ~ ip. Substi- 
tution of these relations into eq. (|5.7[) yields the relation 

s(<p) ~ [^(^l-'ElloO + Wi+dvMv) ~ <p which 

is in agreement with the estimation (|5.10|) [taking into 
account that at ip 3> 1 in eq. (|5.6p y>coth(2<^) « y?]. 



c. Interpolation formula. The limiting results obtained 
above for y> -C 1 and can be combined into a simple 

algebraic interpolation formula 



?n(v) 



2.p 



sinh(2^) 6 + 2ip 2 



(5-11) 



which reproduces function Q n (<p), numerically evaluated 
using eqs. (|5.6p - (|5.8p . with accuracy ~ 3% (see Fig. 2a). 

B. Wide- well shape 

Another form of the well shape, in which analysis of 
the escape rate w e can be made analytically, corresponds 
to the small distance of closest approach, or large a, for 
which d ~ r& <C a. In this case in a fairly wide region of 
relatively strong force / < 1/d, 1/r^ the escape kinetics is 
fairly accurately described by the two-state model (|4.1bp . 

It is important to note that the inequality fr^ < 1 en- 
sures quite high accuracy of the approximation neglect- 
ing the effect of force on the well shape in the region near 
the bottom. In this approximation, the quasicquilibrium 
population distribution within the well is isotropic: 



I*, 



\Yo) 



(5.12) 



This, in turn, means that the partition function Z w is 
independent of the angle 9 and is given by eq. (|3.6p . 

The effect of force, however, manifests itself in the 
anisotropy of the activation energy of escaping u a (9) : 

u a {9) ~ u f (0,a) w u b + 2tpcos9, (<p = fa/2), (5.13) 

which, in turn, leads to the anisotropy of the detailed 
balance relation, i.e. the anisotropy of the equilibrium 
constant 



K e {9) = K° e e- 2vcosl 



(5.14) 



where K® is the isotropic equilibrium constant in the 
absence of external force given by eq. (|5.13p . 

Formula (|5.14p calls for some additional comments es- 
pecially concerning its applicability. The fact is that 
the value of K e (9) at each particular 9 is determined 
assuming local quasiequilibrium of the population out- 
side and inside the well in the region close to r — a at 
this 9. In general, it is difficult to justify the existence of 
the quasiequilibrium in the considered limit, unlike the 
limit of narrow well discussed above. This is because for 
ip = fa/2 < 1 the time of passing over the escaping bar- 
rier width Sh ~ min {a, (the width of the region of 
transition from the inner part of the well to the outer 
one), Tb ~ fil/D is comparable with the time of reorien- 
tation r c ~ a 2 /D. It is worth noting, however, that the 
accuracy of quasiequilibrium assumption becomes better 
with increasing / since the for ip — fa/2 ^> 1 the width 
6b <C a and, correspondingly, Tb <C r c . 

The above-mentioned arguments lead us to the con- 
clusion that in the considered limit of small radius of the 
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well bottom the two-state model with ^-dependent equi- 
librium constant K e (9) gives quite reasonable interpola- 
tion formula for the kinetics of the escaping process and, 
in particular, for the escape rate w e , which correctly de- 
scribes both the limit of weak and strong external force. 
Further analysis (see below) will confirm this statement. 

Formula for the escaping kinetics can straightforwardly 
be derived with the use of general formulas (|4. 39[) , (14.401) , 
and some results obtained above in the limit of narrow 
potential well. The fact is that, in the mathematical 
form, the average of any operator multiplied by angu- 
lar dependent equilibrium constant (|5.14[) [of type of eq. 
(|4.41[) ] over the isotropic equilibrium state (|5 . 12[) is simi- 
lar to the average over the equilibrium distribution (15.4|l . 
except for the partition function Zo(tp) [eq. (|5.5|) ]. which 
should be replaced by Zo(ip — > 0) = 2 corresponding to 
the isotropic distribution. These simple algebraic manip- 
ulations result in the following expression for w e (tp) 

w e (tp)/w° ew (tp) = Q w {ip) = ±Z (<p)e- 2v Q n (<p) (5.15) 

where Z§(tp) and Qniv) are determined in eqs. (|5.5p and 
(|5.6|) . respectively, and 

w° e J<p) = wy (5.16) 

is the escape rate in the absence of the external force 
but with the activation energy u* a , corresponding to the 
orientation 9 = tt (most favorable for escaping): 

u* a = u a (9 = tt) = Ub - 2tp. (5-17) 

As in the case of narrow potential well the dependence 
of w e on the force / is expressed in terms of that on the 
only parameter tp. The characteristic function Q w {tp), 
which determines the pre-exponential factor in the acti- 
vation type dependence of w e {tp), is displayed in Fig. 2b. 
The numerical results show that Q w (tp) monotonically 
decreases (with increasing ip) from Q w = 1 at ip = to 
Q w = 1/2 at ip — > oo. This behavior is markedly differ- 
ent from that of Q n {tp) although the specific features of 
<3uj (^-dependence are essentially determined by those of 
Q n (tp). Some of features of the function Q w (tp), for ex- 
ample saturation at ip — > oo, looking unexpected at first 
sight, can be understood by simple analysis (see below). 

a. Weak force limit. In the weak force limit the behav- 
ior of Q w (ip) at ip = fa/2 -C 1 differs form that obtained 
above for Q n (tp) (i.e. for narrow potential well): Q w (ip) 
decreases with increasing tp, so that at ip <C 1 

QM = Q w (tp<^l)^l-tp. (5-18) 

Such a behavior of Q w (<p) results from using the in- 
dependent normalizing rate bj^ ~ e 2<p (instead of W® = 
Wg) in the definition of Qwi^p)- 

b. Strong force limit. In the opposite limit tp = 
fa/2 ^> 1 the force strongly affects the average escape 
rate w e , first of all, because it significantly changes the ac- 
tivation energy of the rate w e . As for Q w (tp), which char- 
acterizes the pre-exponential factor of the corresponding 



Arrenius-type expression for w e , at tp 3> 1 it monotoni- 
cally decreases approaching the asymptotic value 1/2. 

The obtained Q w ((/^-independence at tp — > oo can eas- 
ily be understood by taking into account that, according 
to formula (|5.14[) . in the case of wide well the strong 
external force causes significant localization of the flux 
of escaping particles in a small region of orientations 
59 = tt — 8 < 1/^/tp <C 1. The escape rate is deter- 
mined by the total flux J e through this region of size 
s e ~ {50) 2 ~ tp^ 1 . In the strong force limit tp ^> 1 
the flux J e ~ ip, as it follows from eq. (|5.10[) . so that 
Qwif) ~ Se Je ~ const. The exact estimation can be ob- 
tained just by substitution of the corresponding limiting 
expression (|5.10|) into eq. (I5.15|) : 

Q s w (tp) = IZoMe-^MUr = \ (5.19) 

Similarly to the narrow well limit, in the case of wide 
well for large tp the escape rate is determined by the 
quasi-one-dimensional flux of escaping particles. The 
mechanism of formation of the one-dimensional flux is, 
however, somewhat different in both cases: for narrow 
wells the transition to the one dimensional regime re- 
sults from high localization of the well population in the 
small region at 9 ~ tt, while for wide wells this transition 
is caused by strong localization of favorable transition 
rates in this region. 

c. Interpolation formula. A simple interpolation ex- 
pression for Q w (ip) can be derived, for example, with the 
use of similar formula for Q n {tp) presented in eq. (|5.11[) : 

QM « Qi(<P) - IZ^e'^Q^). (5.20) 

Quite satisfactory accuracy (~ 3%) of this formula is 
demonstrated in Fig. 2b. 

VI. SUMMARY AND CONCLUDING 
REMARKS 

This work concerns detailed theoretical study of the 
effect of the external force / = F/{kBT) on the kinetics 
of depopulation of a deep potential well. The well is as- 
sumed to be isotropic and highly localized (short range) . 
Though, detailed analysis show a 16 i 17 that fairly deep po- 
tential well resulted, for example, from the Coulomb in- 
teraction (at large r) can also be considered as highly 
localized in some conditions. 

Fairly simple matrix expression for the kinetics of the 
well depopulation is obtained, which predicts strong ef- 
fect of external force on the kinetics. Moreover, in gen- 
eral, the kinetics is predicted to be non-exponential. 

In our work we have concentrated on the analysis in 
the most physically reasonable limit of fast orientational 
relaxation of the PDF in the well. In this limit the ana- 
lytical expression for the depopulation kinetics is derived 
which predicts the kinetics to be close to the exponential 
and the total depopulation rate to be a sum of reaction 
and escape rates. In our work we have mainly studied 



the specific features of the escape rate w e whose value 
appears to significantly depend on shape of the well. Sim- 
ple analytical expressions for w e (f) are obtained for two 
limiting types of wells: narrow wells of type of well layer 
at a distance of closest approach d (for which a — d <C a) 
and deep wells with small distance d -C a. 

In the case of narrow well the effect of the force on the 
escape rate is fairly strong but shows itself only in the 
preexponential factor of the Arrenius-type dependence 
of the rate, i.e. no strong effect on the activation energy 
is predicted. On the contrary, in the case of wide well 
(d <C a) the force affects not only preexponential factor 
but the activation energy as well. 

It is worth noting that the effect of an external force on 
the diffusion-assisted processes in the presence of interac- 
tion between particles are studied in a number of works 
(see, for example, refs. [1] and [7]). Especially compre- 
hensively the force effect (electric field effect) is analyzed 
in the case of ion pair recombination reaction, i.e. in the 
case of the Coulomb interaction between particles. 

Unfortunately it is practically impossible to compare 
the results our analysis with those obtained earlier, since 
the earlier works mainly concerned with processes in po- 
tentials without well at short distances, the reactivity is 
usually assumed to be high. In particular, in the case of 
ion pair recombination processes the recombination ki- 
netics is considered to be determined by diffusive motion 
in the pure Coulomb potential^^ It is, nevertheless, in- 
teresting to note that in the small field limit fa <C 1 the 
force effect on the probability P e (f) of escape from the 
Coulomb potential, found in ref. [4], is independent of 
the initial distance between ions and is represented in the 
form P e (f) ~ Pe(f = 0)(1 + /a/2), which is in apparent 
agreement with the field dependence of the escape rate 
obtained in our works [see eqs. (|5.9[) and (|5.18[) ]. 

Noteworthy is also that in the strong force limit the 
escaping process becomes nearly one-dimensional in both 
cases of well shape considered. In this limit the escaping 
rate is determined by the one-dimensional flux from the 
small region of favorite orientations corresponding to 9 ~ 
7r. This fact allows one to easily improve the considered 
two-state model, in which the effect of the force on the 
location of top of the barrier (assumed to be at r = a) 
is neglected. Moreover, in the strong force limit one can 
also take into account the smoothness of the shape of 
the realistic barrier near the top, which in the two-state 
model is actually assumed to be of cusp shape. 

Concerning possible applications of formulas obtained, 
note that the most convenient for experimental analysis 
is not the force dependent inverse mean lifetime wo((p) 
[see eq. (|4.40[) ] . but the difference wo(ip) — wq(0) = 
w e ((f) — w®, which is independent of the rate w r of reac- 
tion within the well (assumed to be independent of ip). 
The corresponding dimensionless parameters 

SQv((p) = [wo((p)-w (0)]/wl v ((p) ) (v = n,w), (6.1) 

are directly related to Q v ((f): 

6Q n (<p) = Qn(if) ~ 1, 5Q w (<p)=Q w (<p)-e- 2 *. (6.2) 



10 




0,04 ■ 1 ■ 1 ■ 1 . 1 . 1 ■ 

1 2 3 4 5 <p 



FIG. 3: The dependence of the dimensionless difference of 
rates 5Q m (<p) [see eqs. (|6.1[) and (|6.2[l ] calculated with the ex- 
act formula (full lines) and the interpolation expression (|5.20p 
(triangles) . 



The behavior of 8Q n (ip) is, clearly, similar to that of 
Qni}p) except for evident displacement along the ordinate 
axis. As for <5(3™(y)-dependence, shown in Fig. 3, its 
form is essentially different from that of Q w (</?): at <p — > 
the function 5Q w (tp) rs cp is similar to SQ n (tp), while 
SQ w (ip — > oo) = 1/2. Moreover 5Q w (ip) has a maximum 
(though not very pronounced) at ip = <p m ~ 2.0. 

It is of great interest to compare these theoretical pre- 
dictions with experimental results of type of those given 
in refs. [7-14] but in the presence of electric field. 

Concluding our brief discussion we would like to note 
that in this work we restricted ourselves to the analy- 
sis of the most realistic limit of fast orientational relax- 
ation within the well. In reality, however, with the use 
of general formula (|4. 16|) one can also describe the man- 
ifestation of finitcness of the orientation relaxation time. 
The case, in which the effect of finiteness is largest, of 
course, corresponds to L c = 0, i.e. the absence of ori- 
entational relaxation. In this case the angular depen- 
dence of the equilibrium rate [K e {9)], evidently, results 
in the highly non-exponential well depopulation kinetics 
no(£), which can be approximated by the sum of expo- 
nentially decreasing (monomolecular) contributions with 
^-dependent rates, coming from different orientations. 
With the use of obtained formulas there will be no diffi- 
culties to analyze this case as well, when needed. 
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